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Bound states near a moving charge in a quantum plasma 
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It is investigated how the shielding of a moving point charge in a one-component fully degenerate 
fermion plasma affects the bound states near the charge at velocities smaller than the Fermi one. 
The shielding is accounted for by using the Lindhard dielectric function, and the resulting potential 
is substituted into the Schrodinger equation in order to obtain the energy levels. Their number 
and values are shown to be primarily determined by the value of the charge and the quantum 
plasma coupling parameter, while the main effect of the motion is to split certain energy levels. 
This provides a link between quantum plasma theory and possible measurements of spectra of ions 
passing through solids. 



I. INTRODUCTION 

Quantum plasmas have recently received a rapidly 
growing interest [H-Ql- The motivation is related, in 
particular, to nanostructured metallic and semiconduc- 
tor materials 0, Q , laser-plasma systems [1, , astro- 
physics [ll|, and ultracold plasmas Quantum plas- 
mas can be described by various models including the hy- 
drodynamic one with the Bohm potential 0, fisj ]. kinetic 
Wigner-Poisson model (a quantum analogue of the clas- 
sical Vlasov-Poisson system) 0, HI] , quantum relativistic 
plasmadynamics ^ , and other models • Based on 
these and similar approaches, recently there have been 
studies of quantum effects on various plasma phenomena 
including two-stream instability [l5|, ion-acoustic waves 
[l6l [iTj . spin effects on the plasma dispersion [3, fisf - 
whistlers [23|, Bernstein modes [2l|, Zakharov dynamics 
[2^ , and thermodynamic properties [H, [23| ■ 

As a part of research in this field, there have been 
numerous works to link theoretical models to practice. 
For instance, already in 1956, Watanabe showed that 
his measurements of energy losses and scattering angles 
of electrons passing through a thin metallic film are in 
a remarkable agreement with the quantum plasma dis- 
persion relation that can be derived from the Wigner- 
Poisson model [25). As regards recent works, Manfredi 
and Hervieux performed quantum plasma simulations 
that reproduce low frequency nonlinear oscillations re- 
vealed in transient refiection experiments on thin films 
[26j . Furthermore, Marklund et al. demonstrated that 
the quantum collisionless damping of surface plasmons 
poses a fundamental size limit for plasmonic devices f27| . 

Another link between theory and practice is related to 
spectra of ions moving in quantum plasmas. Indeed, a 
plasma shields the Coulomb potential of a moving charge 
and thus affects the bound electron states near it. This 
was demonstrated experimentally in the regime where 
the charge velocity exceeds the Fermi one HI]. In 
this regime, the effect has been studied in a number of 
theoretical papers [29l - [3^ . 

The object of this Letter is to investigate this effect 



for charge velocities smaller than the Fermi one. In this 
regime, as opposed to that where the charge velocity ex- 
ceeds the Fermi one, two factors may be particularly sig- 
nificant. The first one is the spread of velocities in a fully 
degenerate distribution, as the plasma can no longer be 
considered as cold or almost cold. The second one is the 
quantum recoil/tunneling. Its role in the shielding of a 
moving charge was studied in detail in Ref. [ssj . 



II. MODEL 

We consider a point charge +Ze representing an ion 
and moving through a plasma at constant velocity v. 
Here, e is the elementary charge. The plasma is assumed 
to be a fully degenerate electron gas immersed in a neu- 
tralizing background. We adopt the frame of the moving 
ion and use it as our origin of coordinates. 

In describing the quantum-mechanical dynamics of a 
single electron near the ion in the presence of the plasma, 
we consider the response of the plasma to be solely due 
the ion; this is justified when the plasma responds much 
more strongly to the ion than the electron, as will be 
the case [3J| for Z 3> 1. Thus, the problem is reduced 
to solving the time- independent Schrodinger equation for 
the one-electron wavefunction V'(r), namely 



2m, 



(1) 



where rUe is the mass of the electron and ^p{v) is the elec- 
trostatic potential generated by the ion in the plasma. 
Note that, due to the rotational symmetry of the poten- 
tial about the axis of motion, we can impose the depen- 
dence on the azimuthal angle in spherical or cylindrical 
coordinates to be of the form oc e™"^ for some integer 
value of the magnetic quantum number to. The binding 
energies will be identical for positive and negative to. 
The potential ip{v) can be calculated in the linear ap- 
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proximation by: 
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where k) is the dielectric response function of the 
plasma and eo is the permittivity of free space. Here we 
use Lindhard's dielectric function [36l,f37!| for the response 
of fully degenerate electrons, which reads 



D{LU,k) = 1 + 

3ujp F{uj + iiy + a, k) — F{ijJ + iiy ~ a, k) 



where 



a = , 

F(f^,fc)^^ + i^^4l^lnf^ 
^ ' 2 AkvF \Q-kvF 



(3a) 

(3b) 
(3c) 
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the electron Fermi velocity, 

= -y/ nee^/(rneeo) is the electron plasma frequency, 
and 7ie is the electron number density; v represents an 
infinitesimal positive number. The Lindhard response 
function can be derived from the Wigner-Poisson system 
and includes both the quantum tunneling and the degen- 
eracy. The potential resulting from the Lindhard dielec- 
tric function has been studied in Refs. 31, 33, 35]. The 
applicability of our model is analyzed in the discussion 
section. 

We have solved Eq. ^ numerically by finite-difference 
discretization, with the resulting sparse matrix eigen- 
value problem solved using the ARPACK library |38| . 
The potential (^(r) is found by numerical integration of 
Eq. ([2]) . The bound states in this problem are determined 
(up to scaling) by the following dimensionless parame- 
ters: Z, the atomic number of the ion; the ratio v /vf] 
and H — Hujp/ {meVp), which characterizes the strength 
of the quantum recoil and is also related to the plasma 
coupling parameter. 



III. RESULTS 

Note that as _ff ^ 0, the quantum tunneling be- 
comes unimportant to the screening so that the non- 
dimensionalized potential [i.e. with position normal- 
ized to A = vp/i^p, and the potential normalized to 
Ze/{47reoX)] tends towards the "semiclassical" form. The 
"semiclassical" form implies that it can be obtained from 
the classical Vlasov-Poisson system with the fully degen- 
erate velocity distribution. In this limit, the potential 
surrounding a stationary test charge tends to the Debye- 
Hiickel one, with the resulting screening length being 
equal to the Thomas- Fermi length A/\/3. 
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FIG. 1: The number of distinct energy levels (shown by nu- 
merals 0, 1, . . ., 7) as a function of v/vf (horizontal axis) 
and ZH (vertical axis), in the limit -ff — >■ 0, ZH held con- 
stant. The lines indicate the boundaries between the areas 
corresponding to different numbers of distinct energy lev- 
els. The line style reflects the quantum numbers (with the 
I quantum number determined by tracing back to t; = 0) 
of the bound state that appears when the line is crossed: 
solid for Z = 0, m = 0; dashed for I — l,m = 0; dotted 
for I — l,m — ±1. 



Once the non-dimensionalized potential has been spec- 
ified, the bound states depend (up to scaling) only on 
ZH. The physical meaning of the product ZH is that 
it is equal (up to a constant coefficient of the order of 
unity) to X/{ao/Z), where = AneolT? / (niee^) is the 
Bohr radius. Note that uq/Z is the characteristic scale 
of the ground state of the ion in vacuum. Therefore it 
is informative to consider the limit — >■ 0, ZH held 
constant. 

Fig. [1] shows the number of energy levels in the limit 
H —i' 0, ZH held constant. As ZH increases, the ratio 
of the screening length to the dimensions of individual 
bound states becomes larger, and thus more bound states 
appear and the binding energy of levels increases, tending 
to the bare Coulomb value \E\ = 0.5Z/n^ a.u. (where n 
is the principal quantum number) as ZH — ^ 00. The 
splitting of energy levels occurs because the isotropy of 
the potential is broken. Fig. [2] shows how a state which 
at ?; = is isotropic is distorted by the changes to the 
potential. 

In general, the effect of the motion on the potential is 
greater screening (and the development of a repulsive po- 
tential) behind the moving ion (and to a lesser extent in 
the direction perpendicular to the motion), and decreased 
screening in front of the ion [35] . Thus, since states with 
different angular quantum numbers (i.e. I and m) have 
different angular dependencies, they are affected by the 
anisotropy of the potential in different ways, as can be 
seen in Fig.[TJ In particular, some of the boundary curves 
in Fig. [1] cross, with the result that for particular values 
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FIG. 2: The potential generated by the ion, at n = 0.9vf in 
the hmit iif — >■ 0; and one of the bound states in this potential, 
with ZH held constant and equal to 6. The ion is located at 
the origin and moving to the right; the coordinates are in 
units of A = vp/ijJp- (a) The ratio of the potential if{r) to the 
unscreened Coulomb potential Ze/ (Aneor). The contours are 
equally spaced, and the thicker contour corresponds to the 
boundary between positive and negative potential, (b) The 
wavefunction ip{r) for a particular bound state, which can be 
continuously traced back to the bound state near a stationary 
test charge with quantum numbers n = 3, / = m = 0. 

of ZH^ there can be narrow ranges of v/vp in which a 
bound energy level disappears, or an additional level ap- 
pears. However, as states become more tightly bound, 
their spatial extent becomes much smaller than A, and 
they cease to be greatly affected by the motion. 

Interestingly, the minimal value of ZH at which at 
least one bound state exists is almost independent on 
v/vp and is about 0.6. At the same time, the value of 
ZH at which the second level appears shows a more no- 
ticeable dependence on velocity. This value decreases 
from about 2.3 at zero velocity to about 1.8 at the Fermi 
velocity and thus changes by about 20%. 

In the fully quantum case of finite H (but less than 



FIG. 3: The binding energies (divided by the atomic number 
Z) as a function of velocity for different values of H, but all 
with ZH « 3.82. Different lines of the same style correspond 
to different bound states for the same parameters; all possible 
bound states are shown for given parameters. 

unity), the quantum tunneling does not introduce quali- 
tative changes to the above, though it can make notice- 
able quantitative changes. This can be seen in Fig. [31 It 
shows that for a fixed ZH the binding energies normal- 
ized by Z increase with H and tend to their semiclassical 
values at — !■ 0. Thus, the main effect of the quantum 
tunneling on Fig. [T] would be to displace the boundary 
curves down. 



IV. DISCUSSION 

The use of our model is limited because it does not 
include correlations of plasma electrons, nonlinear and 
relativistic plasma effects as well as the plasma response 
to the bound electron. Let us discuss the corresponding 
applicability limits in terms of the parameters H and Z . 

The neglect of correlations of plasma electrons is jus- 
tified when H <^ 1. Indeed, the use of the collision- 
less mean-field approximation requires that the quantum 
plasma coupling parameter is small. The latter is de- 
fined as r = [e^n^/^/(47reo)]/[TOei;^/2] which is the ratio 
of the characteristic potential energy of interaction be- 
tween neighboring electrons to the Fermi energy. The 
parameter F can be expressed as F = i/^(37r^)^/^/(27r), 
so we get the requirement H ^ 1. In the limit — >■ 0, 
however, the quantum tunneling disappears [isj . Thus, 
the fact that our model includes the quantum tunneling 
but does not include particle correlations makes it incon- 
sistent in a certain sense. Nevertheless, it is widely used 
in the literature, i.e. at small H the quantum tunneling is 
believed to be more important than particle correlations. 

The use of the linear response formalism [Eq. ^] re- 
quires H < This limitation can be derived from 
the condition that the potential energy of the interaction 
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of a plasma electron with the test charge at the length A 
is much smaller than the Fermi energy. 

The relativistic plasma effects are negligible when 
H ^ 0.06. This condition follows from the require- 
ment that the Fermi velocity is much smaller than the 
speed of light. The latter requirement can be written as 
H > 2^a/(37r) w 0.06, where a = e^/{Aneohc) w 1/137 
is the fine structure constant. 

The plasma response to the bound electron should be 
unimportant when Z 3> 1. In this case the ion should 
generate a significantly stronger plasma response due to 
its larger charge. The requirement Z 3> 1 is also evident 
from the results of Ref. [SJ] . 

The four above requirements can be easily satisfied at 
small H and large Z. Note that the relativistic restric- 
tion is formally not met in the limit 77 — 0, ZH held 
constant. This means that taking this limit within our 
model gives results that apply at not extremely small, 
though sufficiently small, values of H. 

It follows that metals are at the edge of applicability 
of our model. For instance, for aluminum the parameter 
H is about 0.7, as can be deduced from its Fermi energy, 
11.7 eV. This value of H suggests that nonideal plasma 
effects d, II] can provide a non-negligible contribution. 
Nevertheless, the Lindhard dielectric function can apply 
quite well in metals, as evidenced by the experiment of 
Watanabe [H]. 

V. CONCLUSION 

We have studied how the shielding a moving charge 
in a fully degenerate electron gas affects the bound elec- 



tron states near the charge for velocities smaller than the 
Fermi velocity. We accounted for the shielding by using 
the Lindhard dielectric function. 

The main result is that the number of distinct en- 
ergy levels is primarily determined by the parameter ZH, 
while the main effect of the motion is to split certain en- 
ergy levels. In the semiclassical limit, i.e. j? — >■ 0, ZH 
held constant, the minimal value of ZH at which at least 
one bound state exists is almost independent of veloc- 
ity and is about 0.6, while the value of ZH at which 
the second level appears varies from 2.3 to 1.8 with ve- 
locity. In the fully quantum case of finite H (but less 
than unity), the quantum tunneling does not introduce 
qualitative changes, but it can introduce noticeable quan- 
titative changes. They are that the binding energies nor- 
malized by Z increase with H, for a fixed ZH . 

This provides a link between quantum plasma theory 
and possible measurements of spectra of ions passing 
trough solids. For instance, the measured positions of 
spectral lines as well as their splitting could be compared 
to theoretical values. This might allow, in particular, as- 
sessing the role of the effects not included in the linear 
mean-field model, for various conditions. 
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